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1.0  INTRODUCTION 

The  usual  methods  of  calculation  of  wind  tunnel  wall  Interference  (as  discussed  in 
Refs.  1 through  10)  are  generally  unsatisfactory  for  those  cases  where  extreme  downwash 
is  required  for  lift  generation,  such  as  for  VSTOL  configurations.  A study  was  conducted 
to  evaluate  the  vortex  lattice  technique  as  a tool  for  the  calculation  of  wind  tunnel  wall 
interference.  Vortex  lattice  simulations  were  formulated  for  closed  wind  tunnel  walls  and 
for  both  perforated  and  slotted  boundary  conditions.  Solutions  were  obtained  for 
selected  basic  tunnel  configurations  and  were  compared  with  results  from  other  available 
theoretical  techniques.  In  addition,  solutions  for  a more  complex  wind  tunnel/lifting 
model  combination  were  computed  to  demonstrate  the  ability  to  calculate  interference 
for  configurations  for  which  no  exact  analytic  solution  is  feasible. 

2.0  VORTEX  LATTICE  FORMULATION 

The  vortex  lattice  technique  is  a numerical  method  for  solving  the  three-dimensional 
Laplace  equation.  A computer  program  recently  developed  at  AEDC  is  capable  of 
representing  complex  arbitrary  shapes  such  as  the  one  illustrated  in  Fig.  1,  from  Ref.  11. 
The  program  can  calculate  velocities  and  pressure  coefficients,  both  on  the  body  surface 
and  in  the  flow  field.  Streamlines  and  contour  plots  or  parameters  such  as  flow 
angularity,  pressure  coefficient,  etc.  are  also  available. 

The  program  was  adapted  to  compute  the  classical  lift  interference  factor,  6,  defined 
in  Ref.  1 as 

s.  (c/scL)  f/v.)  (1) 

Results  were  found  to  correlate  well  with  solutions  for  several  basic  closed  wall 
tunnel  configurations  (Refs.  1 and  4).  For  example,  the  distribution  of  interference  along 
the  test  section  centerline  for  a single  horseshoe  vortex  is  presented  in  Fig.  2 for  four 
tunnel  cross  sections.  Included  is  the  distribution  of  interference  computed  using  the 
vortex  lattice  representation  of  the  NASA/Ames  40-  by  80-ft  flat  oval  wind  tunnel  shown 
in  Fig.  3.  A modified  interference  factor,  6',  has  been  used  to  illustrate  the  effect  of 
cross-sectional  area  as  well  as  shape.  The  factor,  defined  as 

S = SC  / C (2) 

relates  the  interference  of  each  tunnel  to  that  of  a square  test  section  having  a 
cross-sectional  area,  C'. 
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b.  Vortex  lattice  representation 
Figure  1.  Typical  application  of  the  vortex  lattice  technique. 
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Figure  2.  Theoretical  wall  interference  for  several  closed  tunnel 
cross  sections. 


Figure  3.  Vortex  lattice  representation  of  the  NASA/Ames  flat  oval  tunnel. 
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Replacement  of  the  closed  wall  with  the  appropriate  boundary  condition  would 
permit  calculations  similar  to  those  presented  above  for  wind  tunnels  with  slotted  walls. 
The  development  of  the  necessary  boundary  condition  is  described  in  the  following 
sections.  Details  of  the  numerical  method  used  to  obtain  solutions  are  documented  in 
Ref.  11. 

2.1  SOLID  SURFACES 

The  vortex  lattice  method  uses  an  arrangement  of  vortex  singularities  to  define  solid 
bodies  and  other  geometrical  forms  such  as  those  illustrated  in  Figs.  1 and  3.  For  each 
vortex  element  of  strength,  H,  one  boundary  condition  may  be  specified  at  a point 
(control  point)  in  the  field.  Generally  a control  point  is  located  at  the  centroid  of  each 
vortex  ring  or  horseshoe.  To  define  a ,solid  body  surface  the  boundary  condition  imposed 
at  each  control  point  restricts  the  local  induced  velocity  so  that  it  is  parallel  to  the 
surface  at  that  point. 

The  numerical  method  involves  describing  the  vortex  network  with  a system  of 
linear  equations  which  are  solved  for  the  vortex  strengths. 

For  a general  system  of  vortex  elements,  the  velocity  at  the  ith  control  point  is, 
from  Ref.  1 1 , 


V : = V . + V 


Vi 


V 


N 

z 

i=i 


r\ 


cu 


(3) 


where  V is  the  free-stream  velocity  and  Vvi  is  the  velocity  induced  at  the  control  point 
by  the  vortex  system.  The  influence  coefficient,  Cy,  is  the  velocity  induced  at  the  ith 
control  point  by  the  jth  singularity  when  its  strength  is  unity. 

The  most  common  form  of  boundary  condition  used  in  the  vortex  lattice  method  is 


To  satisfy  the  expression  at  the  ith  control  point  the  component  of  the  local  velocity,  V), 

A 

in  the  direction  parallel  to  the  unit  normal  vector,  bj,  must  have  a magnitude,  Gj.  For  a 
solid  boundary  the  velocity  normal  to  the  surface  is  zero.  Combining  Eqs.  (3)  and  (4) 
yields,  for  all  control  points, 


V, 


A 

= V. 


r.  = 

j 


= 1,2.. .N 


(5) 
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Rearranging, 


2(c1J.b,)rj  . -v_.  (6) 

The  set  of  N linear  simultaneous  homogeneous  equations  can  be  solved  to  yield  the 
singularity  strengths  (rj,  i = 1,2  . . . N). 

The  velocity  at  any  point,  F,  in  the  flow  field  can  then  be  calculated  as  follows. 
Influence  coefficients  for  the  particular  point  are  determined  by  using  the  Biot-Savart 
equation  to  compute  the  velocity  induced  by  the  straight  line  segment  vortices  of  each 
vortex  loop.  Combining  the  influence  coefficients  provides  a matrix,  Q.j,  which  can  be 
substituted  into  Eq.  (3),  yielding 


£ rj?n 

i=l 


The  pressure  coefficient  can  be  computed  at  the  same  point,  using 


2 


(7) 


(8) 


Computation  of  either  velocity  or  pressure  coefficient  in  the  immediate  vicinity  of 
the  vortex  sheet  requires  special  consideration  because  of  the  close  proximity  of  the 
finite  singularities.  Representative  surface  velocities  are  available  only  at  the  control 
points  since  the  surface  boundary  condition  is  not  satisfied  elsewhere. 

In  determining  the  velocity  on  either  side  of  the  surface  at  a control  point,  one  must 
account  for  the  velocity  jump  across  the  discontinuity  sheet.  The  jump  velocity  is  equal 
to  the  cross  product  of  the  vortex  sheet  density  vector,  Kj,  and  the  unit  vector  normal  to 
the  local  surface.  The  sign  convention  is  illustrated  graphically  in  Fig.  4.  For  the  ith 
control  point  the  difference  between  the  outer  and  inner  surface  velocities  can  be  written 


VDi  “ Ki 


A 

x b, 


(9) 


The  unit  vector,  Bj,  is  a special  case  of  the  unit  normal,  bj,  and  is  included  to  define 
the  sign  of  the  velocity  change.  Due  to  the  discrete  nature  of  the  lattice-type  model, 
some  care  must  be  taken  to  evaluate  the  vortex  sheet  density,  Kj.  A technique  similar  to 
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the  one  described  by  Rubbert  in  Ref.  12  was  used  in  this  study.  The  procedure  involves 
dividing  each  individual  vortex  filament  adjacent  to  the  control  point  by  an  appropriate 
length  to  yield  a density.  The  individual  densities  are  averaged  to  provide  an  effective 
local  sheet  density  in  terms  of  the  unknown  singularity  strength,  Ti. 


b.  Application  of  velocity  correction 
Figure  4.  Graphic  representation  of  the  jump 
velocity  correction  technique. 

The  velocity  can  be  determined  on  the  inner  side  of  the  singularity  sheet  by  adding 
one  half  of  the  jump  velocity  from  Eq.  (9)  to  the  velocity  induced  by  the  free  stream 
and  the  vortex  network  from  Eq.  (3).  The  inner  and  outer  surface  velocity  at  the  ith 
control  point,  respectively,  can  be  written  as 


% V 


Di 


(10) 
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and 


Figure  4b  graphically  illustrates  the  vector  addition  of  the  jump  velocity  correction 
to  the  velocity  induced  by  the  vortex  system.  The  free-stream  velocity  has  been  omitted 
for  clarity. 

The  inner  and  outer  surface  pressure  coefficients  can  now  be  calculated  by 
substituting  the  velocities  from  Eqs.  (10)  and  (11),  respectively,  into  Eq.  (8). 

2.2  PERFORATED  WALLS 

Several  theoretical  boundary  condition  formulations  have  been  developed  for  dealing 
with  various  wind  tunnel  wall  configurations.  Of  primary  interest  are  those  discussed  in 
Refs.  1 through  4 for  representing  homogeneous  perforated  and  slotted  boundaries.  The 
interference  results  for  rectangular  tunnels  with  ventilated  upper  and  lower  walls  and 
closed  side  walls  have  been  used  throughout  the  present  study  for  comparison  with 
existing  theoretical  calculations.  In  each  case  a single  horseshoe  vortex  of  unit  strength 
was  utilized  to  represent  the  lifting  model. 

The  perforated  wall  boundary  condition  presented  in  Ref.  2 is  intended  to  represent 
a tunnel  wall  with  a uniform  distribution  of  holes.  The  theory  was  developed  based  upon 
an  assumption  that  such  a wall  behaved  in  the  same  manner  as  a porous  material 
boundary  (i.e.,  that  the  flow  through  the  surface  was  linearly  proportional  to  the  pressure 
difference  across  the  boundary).  The  boundary  condition  can  be  written  as  follows; 

£,ci*  - 0 (12) 

ox  an 

Expressed  in  terms  of  the  velocity  on  the  top  wall  of  a wind  tunnel,  Eq.  (12)  becomes 

u + Cw  = 0 (13) 

The  boundary  condition  states  that  the  perturbation  velocity  components  at  the  wall  are 
related  by  a constant,  C.  The  wall  porosity  parameter,  Q,  is  conventionally  defined  as 

Q = 1/(1  r C)  (14) 
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For  a closed  wall  where  Q = 0 (or  C = »),  Eq.  (13)  reduces  to 

w = 0 (15) 

For  Q = 1 (or  C = 0)  the  boundary  condition  represents  a constant  pressure  or  open  jet 
boundary  when 

u = 0 (16) 

The  form  of  the  perforated  wall  boundary  condition  is  similar  to  that  used  for  the 
solid  surface  [Eq.  (4)]. 

Equation  (13)  states  that  the  velocity  in  a particular  direction  (i  + Ck)  must  be  zero 
along  the  top  wind  tunnel  wall.  As  was  the  case  with  the  solid  boundary,  the  scalar 
velocity,  Gj,  is  set  equal  to  zero.  The  unit  normal  vector,  bj,  in  Eq.  (4)  is  replaced  by  a 
unit  vector,  dj,  in  a direction  determined  by  the  wall  porosity.  For  an  upper  tunnel  wall 
the  new  unit  vector  can  be  written 


I.  = 


C It 


V i + c2  \/i+  c‘ 


(17) 


The  solid  surface  formulation  dealt  only  with  restricting  the  velocity  in  the  direction 
normal  to  the  surface.  With  the  introduction  of  the  other  component  (i.e.,  the  direction 

A 0 

k),  it  becomes  necessary  to  account  for  the  effects  of  local  singularity  sheet  density.  The 
velocity  expression  developed  for  computing  the  velocity  on  the  inner  surface  at  h 
control  point,  Eq.  10,  was  used  since  it  contains  the  required  terms.  The  perforated  wall 
boundary  condition  on  the  inner  surface  of  a vortex  lattice  wind  tunnel  becomes 


is 


di  = 


(18) 


The  boundary  condition  is  applied  to  Eq.  10,  yielding,  for  the  ith  control  point, 


V 


IS; 


A 


A 


K;xBi 


(19) 


Evaluating  the  terms  in  Eq.  (10)  for  each  of  the  control  points  at  which  the 
perforated  boundary  condition  is  to  be  imposed  (i.e.,  i = 1,2,...  Np)  yields  Np  linear 
simultaneous  equations  in  terms  of  the  N unknown  Tj  values.  The  remaning  equations,  if 
N,  may  represent  solid  surface  control  points  on  nonperforated  portions  of  the 
tunnel  wall  or  on  the  test  article. 
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Comparisons  of  centerline  wall  interference  distributions  from  Refs.  1 and  3 with 
data  calculated  using  the  vortex  lattice  model  are  presented  in  Fig.  5 for  several  porosity 
parameter  values.  The  agreement  of  the  vortex  lattice  results  with  solutions  from  available 
theoretical  calculations  is  good. 

2.3  SLOTTED  WALLS 

Early  slotted  wall  models  were  based  upon  discrete  geometric  representations  of 
slots  and  solid  panels  (Ref.  13).  A typical  example,  presented  in  Fig.  6,  has  one  wide  slot 
in  each  comer  of  the  cross  section.  The  slots  are  simulated  by  vortex  elements  which  are 
the  same  width  as  the  slots.  A free  jet  or  constant  pressure  boundary  condition  is 
imposed  at  the  slot  control  points  to  permit  flow  through  the  wall.  The  solid  panels  are 
represented  by  vortex  loops  with  conventional  dosed  wall  boundary  conditions  specified 
at  their  centroids.  The  technique  works  well  for  tunnels  with  small  numbers  of  large 
slots,  as  discussed  in  Ref.  13.  For  wind  tunnels  with  several  narrow  slots  a large  number 
of  vortex  elements  is  required  since  the  singularity  spacing  must  be  kept  in  the  order  of 
the  slot  width.  This  in  turn  places  prohibitive  demands  upon  the  computer  capability. 
Since  most  existing  slotted  tunnels  were  designed  with  narrow  slots,  the  discrete 
simulation  technique  was  abandoned. 


a.  A = 1,  t = 0.05,  solid  side  walls,  top  and  bottom  perforated 
Figure  5.  The  vortex  lattice  solution  and  other  calculated  results  for  various 
perforated  wall  tunnels. 


fc 
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x/jSb 


b.  X = 1,  r = 0.05,  all  walls  perforated 


c.  X = 1,  r = 0.6,  solid  side  walls,  top  and  bottom  perforated 
Figure  5.  Concluded. 
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; 


Figure  6.  Vortex  lattice  representation  of  a wind  tunnel  with  discrete  slots. 


The  alternative  method  involved  application  of  a homogeneous  slotted  wall 
boundary  condition  over  the  entire  slotted  wall.  Such  a boundary  condition  has  been 
used  by  several  investigators  (e.g.,  Refs.  1,  5,  6,  7,  8,  9,  and  10).  The  boundary  condition 
attributed  to  Buseman  in  Ref.  6 can  be  expressed  as 


d<fi  d2$ 
~d*  + K dx  dn 


(20) 


where  x is  the  direction  of  the  tunnel  axis  and  n is  outward,  normal  to  the  wall  inner 
surface.  The  geometric  slot  parameter,  k,  is  related  to  the  physical  tunnel  dimensions 
shown  in  Fig.  7 by  Chen  and  Mears  (Ref.  5)  as 


(21) 


and 


P = 1 / [l  - («A)] 


(22) 
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Figure  7.  Slotted  wall  tunnel  cross-sectional  geometry  and  coordinate  system. 

The  slot  parameter,  P,  is  most  frequently  used  since  its  limiting  values  are  0 and  1 for 
closed  and  open  jet  boundaries,  respectively,  whereas  k varies  from  0 to  infinity. 

Written  in  terms  of  the  upper  wind  tunnel  wall,  the  boundary  condition  becomes 

U 4-  « ^ = 0 (23) 

(?X 

where  9w/dx  is  the  local  gradient  of  the  w-component  of  velocity  in  the  x-direction. 
Rewritten  in  terms  appropriate  for  a discretized  representation  with  a uniform  singularity 
spacing,  Ax,  in  the  x-direction,  the  boundary  condition  at  the  ith  control  point  is 

(w.+i  ”*.-■) 


U-  +■  K 


0 


(24) 
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Rearranging, 


Wj- 


i-1 


) 


= 0 


(25) 


Defining  k*  - k/2Ax  yields 


ui  - K*  (wi+l  ” wi-i)  = 0 (26) 

which  is  the  form  of  Eq.  (13).  A unit  direction  vector  analogous  to  Eq.  (18)  can  be 
formulated  as 


A 

i 


4« 

+ K 


'/l  + K 2 


(27) 


By  the  same  analogy,  a control  vector  expression  is  needed  which  is  similar  to  the  inner 
surface  velocity  vector  in  Eq.  (10)  but  modified  to  contain  the  new  parameters 
introduced  by  the  slotted  wall  boundary  condition.  On  the  basis  of  the  direction  vector 
.already  defined  in  Eq.  (27),  such  an  expression  can  be  written  as 


Vfs.  = vw  + Vfc  + V,VDi  = v£r  j C*  + (*K,  X B;  ) 

1 i=i 

In  terms  of  the  wind  tunnel  upper  wall,  the  vortex  system  induced  velocity, 

-*■  * * * 

VVi  = Uj  i + Vj  j + W; k 

has  been  replaced  with  a new  vector, 


(28) 


(29) 


Vi 


ui  i + vi j + 


(- 


i+1 


(30) 


It  is  not  necessary  to  alter  the  free-stream  or  jump  velocity  terms  since  neither  has  a 

component  in  the  z-direction.  The  boundary  condition  can  be  written  as 

I 


(31) 
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which  yields 


(32) 


Evaluating  Eq.  (32)  at  each  control  point  on  the  slotted  walls  and  combining  the 
resulting  equations  with  those  for  other  types  of  boundaries  yields  the  usual  set  of  N 
simultaneous  linear  equations  in  terms  of  the  unknown  vortex  strengths,  Tj. 

Calculations  were  made  for  a basic  slotted  wind  tunnel  configuration  similar  to  the 
one  previously  used  in  the  perforated  wall  analysis.  The  test  section  was  square  in  cross 
section  with  closed  side  walls  and  slotted  upper  and  lower  walls.  Lift  was  generated  by  a 
single  horseshoe  vortex.  Results  of  the  computations  are  presented  in  Fig.  8 along  with 
analytic  solutions  from  Ref.  3.  As  can  be  seen,  correlation  is  excellent  throughout  the 
slot  parameter  range. 


*//?  b 


Figure  8.  The  vortex  lattice  solution  and  other  calculated  results  for  several 
values  of  the  homogeneous  slot  parameter,  P. 
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3.0  APPLICATION  OF  THE  METHOD 

The  vortex  lattice  method  can  predict  lift  interference  in  tunnels  with  closed, 
perforated,  or  slotted  walls  for  lifting  models  (Figs.  5 and  8).  The  configurations 
presented  thus  far  have  been  limited  to  those  involving  an  elementary  horseshoe  vortex  of 
fixed  strength  to  represent  test  article  lift.  The  method  is  equally  valid  for  a model  with 
distributed  lift.  Since  the  formulation  is  linear,  additional  specified  vortex  elements  of 
fixed  strength  can  be  included  to  represent  the  desired  lift  distribution.  Calculations  can 
readily  be  performed  for  test  article  configurations  so  complex  that  analysis  by  other  ✓ 
means  is  impractical. 

The  vortex  lattice  technique  is  also  capable  of  accounting  for  the  mutual  interaction 
between  the  wind  tunnel  walls  and  the  test  article.  The  requirement  of  the  fixed  strength 
lifting  element  distribution  can  be  relaxed  and  replaced  with  a set  of  boundary  conditions 
and  vortices  of  unspecified  strength.  The  solution,  obtained  by  simultaneously  satisfying 
the  tunnel  wall  and  test  article  boundary  conditions,  contains  influence  of  both  the 
tunnel  and  the  test  article.  The  tunnel  wall  effects  can  be  assessed  directly  through  a 
comparison  of  the  mutual  interference  results  and  the  solution  obtained  in  an  infinite 
free  stream  (i.e.,  an  interference-free  environment).  The  following  sections  demonstrate 
the  use  of  unspecified  strength  vortices  to  represent  the  test  article. 

3.1  DISTRIBUTED  LIFT  MODEL 

A more  realistic  model  than  the  simple  horseshoe  vortex  model  was  selected  to 
investigate  the  compatibility  of  the  ventilated  wall  boundary  conditions  with  a distributed 
lift  model  represented  by  unknown  strength  vortex  elements. 

Dimensions  of  the  wing-tail  configuration  selected  are  shown  in  Fig.  9.  The  vortex 
lattice  representation  of  the  lifting  model  is  shown  in  Fig.  10  installed  in  the  30-  by 
45-in.  vortex  lattice  wind  tunnel  used  in  the  investigation.  The  test  section  represents  a 
low  speed  wind  tunnel  at  AEDC  having  the  same  dimensions. 

The  wing  model  is  divided  into  six  sections  along  the  span.  Each  section  is  made  up 
of  200  vortex  elements,  distributed  chordwise  over  the  airfoil  surface,  with  100  on  the 
upper  surface  and  100  on  the  lower  surface.  Both  the  vortices  and  control  points  are 
distributed  by  a cosine  function  to  provide  closer  spacing  near  the  leading  and  trailing 
edges.  The  unit  normal  vectors  located  at  the  control  points  are  illustrated  in  the  same 
figure  along  with  a detail  of  the  wingtip  vortex  elements.  The  vortex  lattice  model  of  the 
tail  was  essentially  the  same  as  the  wing  except  for  scale. 
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CENTER  ON  TUNNEL  CENTERLINE 
Figure  9.  Dimensions  of  the  lifting  model. 

The  wind  tunnel  is  represented  by  equally  spaced  rows  of  vortex  elements  oriented 
parallel  to  the  tunnel  axis,  with  six  rows  on  both  the  upper  and  lower  walls  and  four  on 
each  side  wall.  All  four  walls  were  modeled  to  permit  selection  of  either  the  closed  or 
slotted  boundary  condition. 

3.2  WING/TAIL  PRESSURE  DISTRIBUTIONS 

The  chordwise  distribution  of  pressure  coefficient  was  calculated  for  the  upper  and 
lower  surfaces  of  the  vortex  lattice  wing  and  tail  representations. 

The  interference-free  pressure  distributions  are  presented  in  Fig.  1 1 for  five 
incidence  angles  to  illustrate  the  effect  of  angle-of-attack  variation.  A comparison  of  the 
a - 2 and  -2  deg  data  indicates  that  the  vertical  asymmetry  of  the  wing/tail 
configuration  had  only  a slight  effect  on  the  tail  pressure  distribution.  The  effect  on  the 
wing  pressures  was  almost  undetectable. 
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Figure  10.  Vortex  lattice  representation  of  a wind  tunnel/lifting  model. 


Three  variations  of  the  Fig.  10  configuration  were  investigated,  namely  the  lifting 
model  in  a free  stream  (with  no  tunnel),  the  lifting  mode)  in  a closed  tunnel,  and  the 
lifting  model  in  a tunnel  with  closed  side  walls  and  slotted  upper  and  lower  walls.  In  all 
three  cases  pressures  were  obtained  at  the  quarter-span  location  of  each  lifting  surface. 
The  resulting  pressure  distributions,  computed  at  a = 6 deg,  are  presented  in  Fig.  12. 
As  expected,  the  effect  of  the  closed  wall  test  section  was  to  increase  the  upper  surface 
velocities  and  to  decrease  those  on  the  lower  surface  of  both  the  wing  and  tail.  Thus,  the 
upper  surface  pressure  coefficients  are  more  negative  for  the  closed  wall  case  than  those 
in  the  free  stream.  Since  the  lower  surface  pressures  become  more  positive  (or  less 
negative  over  part  of  the  tail),  a net  increase  in  the  lift  coefficient  would  result.  The 
overall  effect  on  the  pressure  distributions  caused  by  confinement  of  the  flow  field 
within  the  solid  boundaries  appears  to  be  equivalent  to  that  resulting  from  an  increase  in 
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the  model  angle  of  attack.  Such  an  effect  is  consistent  with  theory  (Refs.  1 and  3), 
which  predicts  a positive  interference  in  the  vicinity  of  a single  concentrated  lift 
disturbance  in  a closed  tunnel.  The  interference  upwash  velocity  distribution  was 
computed  using  the  vortex-lattice  closed  wall  tunnel  and  wing/tail  models.  The  results, 
presented  in  Fig.  13,  show  that  a positive  upwash  is  also  generated  at  both  the  wing  and 
tail  locations  for  the  distributed  lift  case. 


Pressure  distributions  computed  within  the  30-  by  45-in.  wind  tunnel  with  slotted 
upper  and  lower  walls  and  closed  side  walls  are  also  presented  in  Fig.  12.  A slotted  wall 


1 */cW  X/C| 

a.  Wing  b.  Tail 

Figure  11.  Pressure  distributions  for  the  wing/tail  model  computed 
in  a free  stream. 
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parameter  value,  P = 0.4,  was  used  in  the  computations.  The  value,  based  on  Ref.  3,  was 
selected  to  provide  minimum  wall  interference  on  the  wing.  The  excellent  agreement  of 
the  calculations  for  the  wing  in  the  slotted  tunnel  with  the  free-stream  solution  {Fig.  12a) 
indicates  that  the  interference  induced  by  the  walls  was  small  at  the  wing  location.  A 
similar  comparison  of  the  tail  pressure  distributions  (Fig.  12b)  reveals,  however,  that  the 
data  do  not  correlate  as  well.  It  appears  that  the  slotted  wall  provided  too  much  flow 
relief,  resulting  in  an  over-correction  in  the  region  of  the  tail. 


a.  Wing  b.  Tail 

Figure  12.  Pressure  distributions  for  the  wing/tail  model  computed 
in  a free  stream  and  in  a 30-  by  45-in.  tunnel. 
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Figure  13.  Normalized  upwash  velocity  computed  along  the  centerline 
of  a 30'  by  45-in.  tunnel  with  closed  walls. 


4.0  CONCLUDING  REMARKS 

A vortex  lattice  technique  for  simulating  the  effects  of  perforations  and 
longitudinal  slots  in  wind  tunnel  walls  has  been  evaluated.  The  lift  interference  factor,  5, 
was  computed  for  a single  lifting  element  in  several  basic  wind  tunnel  configurations. 
Good  correlation  has  been  shown  between  the  vortex  lattice  results  and  other  theoretical 
calculated  solutions.  The  method  was  also  applied  to  a test  article  having  distributed  lift. 
The  wing/tail  configuration,  represented  by  an  arrangement  of  boundary  conditions  and 
vortices  of  unknown  strength,  is  described.  Computed  wind  tunnel  wall  effects  on  the 
lifting  model  pressure  distribution  appear  reasonable. 
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Coupled  with  the  ability  of  the  vortex  lattice  technique  to  accurately  represent 

conventional  test  articles,  the  ventilated  boundary  condition  formulation  has  provided  a 

powerful  tool  for  predicting  the  influence  of  most  perforated  or  slotted  wind  tunnel 

configurations. 
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NOMENCLATURE 


a Slot  width,  ft 

b Tunnel  half  width,  ft 

B,  Unit  normal  vector  at  the  ith  control  point,  negative  away  from  the  boundary 

on  the  inner  surface 

Unit  normal  vector  at  the  ith  control  point,  unit  control  vector 

C Wind  tunnel  cross-sectional  area,  ft2 : constant  related  to  wall  porosity,  Eq.  ( 1 2) 

C'  Cross-sectional  area  of  a reference  wind  tunnel  where  h - 1 and  b = 1,  ft2 

Cjj  Influence  coefficient  matrix,  velocity  induced  at  the  ith  control  point  by  the 

jth  singularity  when  its  strength  is  unity 

Cl  Lift  coefficient,  Cl  = Lift/(q„  S) 

Cp  Pressure  coefficient.  Eq.  (8) 

c Mean  chord,  ft 

d Slot  spacing,  ft 

dj  Unit  vector  at  the  ith  control  point  with  direction  related  to  wall  porosity,  Eq. 

(17) 

ej  Unit  vector  at  the  ith  control  point  with  direction  related  to  the  slot 

parameter,  Eq.  (27) 
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Gj  Magnitude  of  the  component  of  velocity  at  the  ith  control  point  parallel  to  the 

local  unit  control  vector,  Eq.  (4) 

h Wind  tunnel  half  height,  ft 

ij,k  Unit  velocity  vector  components  related  to  the  x,y,z  coordinates,  respectively 
Kj  Vortex  sheet  density  vector,  Eq.  (9) 

M Mach  number 

N Number  of  independent  linear  equations:  number  of  singularities  of  unknown 

strength 

n Coordinate  normal  to  the  local  boundary,  positive  out,  Eq.  (12) 

P Homogeneous  slot  parameter,  l/[l(k/h)] 

Q Porosity  parameter,  1/(1  + C) 

q_  Free-stream  dynamic  pressure,  psf 

S Wing  area,  ft2 

s Wing  semispan,  ft 

u,v,w  Velocity  component  magnitudes  related  to  the  x,y,z  coordinates,  respectively, 
ft/sec 

V Velocity  vector,  ft/sec 

x Coordinate  along  tunnel  axis,  positive  downstream,  ft 

y Horizontal  coordinate  across  tunnel,  positive  left  looking  downstream,  ft 

z Vertical  coordinate,  positive  up,  ft 

a Angle  of  attack,  positive  up,  deg 

(3  Compressibility  factor,  (1  - 

T Vortex  strength  density,  ft/sec 

Ax  Vortex  lattice  spacing  in  the  x-direction,  ft 
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5 Wall  interference  factor 

S'  Modified  wall  interference  factor,  S(C'/C ) 

K Geometric  slot  parameter,  Eq.  (20) 

k*  Modified  geometric  slot  parameter,  k*  = k/2Ax 

\ Tunnel  height -to-width  ratio,  h/b 

r Wingspan-to-tunnel  width  ratio,  s/b 

(J>  Velocity  potential 

SUBSCRIPTS 

D Discontinuity 

F General  point  in  the  Field 

ij  Indices 

IS  Inner  surface 

OS  Outer  surface 

P Perforated  wall,  pressure 

t tail 

V Vortices 

w Wing 

- Free  stream 

SUPERSCRIPTS 
' Related  to  C' 

* Based  on  k* 
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